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In this short note, we investigate the existence of universal relations between the gravimagnetic 
Love number of irrotational stars and the dimensionless moment of inertia. These Love numbers take 
into account the internal motion of the fluid, while the star is globally irrotational. The goal is to ex¬ 
tend the so-called I-Love-Q relations - providing a strong correlation between the gravitoelectric Love 
number, the dimensionless moment of inertia and the dimensionless rotation-induced quadrupole 
- to the gravitomagnetic sector, where internal motion is taken into account. As a byproduct, we 
present for the first time this new gravitomagnetic Love number for realistic equations of state. 


I. INTRODUCTION 

Love numbers encode the deformation of self- 
gravitating fluids towards external tidal perturbations. 
While tidal interactions have been well understood for 
more than a century in the Newtonian context, [1], their 
relativistic counterpart has only started being under¬ 
stood recently. First studies investigated the adiabatic 
limit of static configurations [lE],i.e. under the assump¬ 
tion that the external tidal field varies slowly. Recently, 
the formalism has been extended to slowly rotating con¬ 
figurations HIS], still in the adiabatic limit. A formalism 
to extend the tidal response beyond the adiabatic limit, 
relying on effective field theory has been introduced in 
the Newtonian limit in and extended to the relativis¬ 
tic case in [7]. 

In the relativistic case, static configurations are af¬ 
fected by two sector of gravitational perturbations: odd 
and even, corresponding to gravitomagnetic and grav¬ 
itoelectric perturbations. The gravitomagnetic sector is 
particularly relevant to the case where rotation is present. 
This is a purely relativistic effect and doesn’t have a New¬ 
tonian analogue: a central body surrounded by an orbit¬ 
ing companion will enventually be spun up, or encounter 
internal motion, precisely because of the gravitomagnetic 
field generated by the companion body. 

Gravitomagnetic tidal interactions were first investi¬ 
gated in the post-Newtonian framework in [8] , and in the 
relativistic regime in HE]. However, it has been recently 
shown that taking into account the internal motion of the 
fluid driven by the gravitomagnetic perturbation leads to 
dramatically different values than the gravitomagnetic 
Love numbers assuming a strict hydrostatic equilibrium 
[H]. It was argued in [5] that this assumption of strict 
equilibrium, i.e. no internal motion, is too restrictive for 
realistic situations. 

The value of the tidal coefficients depends on the equa¬ 
tion of state of the fluid. However, it was shown in 
m that gravitoelectric Love number of neutron stars 
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enjoy universal properties, i.e. they are strongly corre¬ 
lated with the moment of inertia and the rotation induced 
quadrupole. These universal relations, dubbed I-Love-Q 
relations involve dimensionless Love number, moment of 
inertia and quadrupole. This result drew a lot of at¬ 
tention, in particular, the universality of the rotation 
induced multipoles have been extended to higher mul¬ 
tipoles mHia leading to an approximate no hair con¬ 
jecture for neutron stars. The validity of the universal 
relations has been extended to the case of rapid rotation 
in [H US], for frequency corrected Love numbers |16j 
and magnetized neutron stars m- These unilversal rela¬ 
tions have further been extended some class of alternative 
models of gravity unmHHni, with the conclusion that 
they remain undistinguishable from predicted by general 
relativity, in the regime already constrained by observa¬ 
tion, with the exception of the dynamical Chern Simons 
theory [52], though its viability can be questioned in the 
non-perturbative regime [23] . These universal relations 
have been recently used for incorporating tidal effects in 
waveform simulations in the post-Newtonian framework 
in [24] . Importantly, universal relations between vari¬ 
ous types of Love numbers (shape, electric and magnetic) 
have been investigated in [25] . In particular, it was found 
that the quadrupole gravitoelectric and gravitomagnetic 
Love numbers (with strict hydrostatic equilibrium) are 
correlated with deviation of the order of 10% for realistic 
equations of state. 


In this paper, we reconsider the universal properties 
of the quadrupolar gravitomagnetic Love number, under 
the assumption of internal motion in the star, along the 
line of [9]. In Sec. |Hj we review the model and equations 
for describing slowly rotating neutron stars and gravito¬ 
magnetic perturbation of neutron stars. In Sec. we 
present our results for the gravitomagnetic Love number 
universality. In particular, we show that a) the gravito¬ 
magnetic Love numbers of irrotational bodies are corre¬ 
lated with the rescaled moment of inertia and b) that the 
universality is enhanced for irrotational bodies, in com¬ 
parison to strict hydrostatic equilibrium. We summarize 
our results in Sec. m 
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II. SLOWLY ROTATING STARS AND 
GRAVITOMAGNETIC LOVE NUMBER 

We use the following line element 

ds^ = -b{r)de + -^ + + sin2 6{d^ (1) 

/W 

+ er{ui{r) - Woo)) + ^■m{ht{r)dt + hr{r)dr)Yidx\ 

where and tm are respectively the small expansion 
parameters for slow rotation and tidal deformation, Woo 
is the pulsar frequency measured by a distant observer, 
Yi is the odd parity vector harmonic, given by Yi = 
{--^d^,sh\Od 0 )Yim{O,ip), where Y/m is the standard 
scalar spherical harmonic, and cc® = (0, Lp). 

We model the neutron star fluid by a perfect fluid 
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FIG. 1. The gravitomagnetic Love numbers for irrotational 
stars (negative) and strict equilibrium (positive), for realis¬ 
tic tabulated equations of state, and three polytropes, as a 
function of the compactness. 


Tab = [pY P)UaUb +Pgab, (2) 


where g is the metric, p is the energy density, P is the 
pressure and u is the four velocity, given by 

= ^(1,0,0, er(w - Woo))® (3) 

VO 


up to order e], and 

We expand the Einstein equations to first order in 
with Cm = 0 for slowly rotating neutron stars, and to 
first order in with = 0 for gravitomagnetic tidal 
deformations. 

The slow rotation equations are given by |26j : 


, „ , 2bUTrGr^P + M) , , h' 

M' = AttGt^p, b' = ^ P' = -(P + p)-, 


n 

LU = 


rir — 2M) 

Auj'{irGr'^ {P + p) — f) 167rGw(P-b p) 


rf 


f 


(4) 


We follow [5] for computing the gravitomagnetic per¬ 
turbations equations leading to 


r^h'' + Vh[ + V{r)ht = 0, (5) 


where 

P ^ 47rr3(P(r) -bp(r)) 

fir) 

^ -2f(r) + 87rgr^(P(r) + p(r)) - (I + 2)(l - 1) 

fir) 

where g = —1 in imiiH]. Recently, Landry and Poisson 
have shown that g = 1 corresponds to the perturbation 
of a non rotating fluid with internal motion [^, named 
irrotational fluid. Here we consider both cases. 

In the following, we will work in geometric units, 
G = c = 1. The vacuum solution to the slow rotation 
equations is given by 


where I is the moment of inertia and Mg the mass of the 
star. 

The vacuum solution of gravitomagnetic perturbation 
equation is expressed in terms of hypergeometric func¬ 
tions: 


hi = 




3(^-1) 


A 4 - 

l ‘ r 


( 8 ) 


where 


. f 2M\ 

A4 = 2F1 I-I + 1,-1- 2; -21; — j , 

B 4 = 2 P 1 — 1, Z -|- 2; 21 -|- 2;-^ . (9) 

We define the following dimensionless quantities 


where G is the compactness, and Rg is the radius of the 
star. Note that [3] defines the compactness as 2Mg/Rs. 

Our goal is to devise a correlation between I and the 
gravitomagnetic Love number, in particular, we will focus 
on the case 1 = 2. We use tabulated equations of state 
for describing the interior of the Neutron Star. 

Before doing so, we present for the first time the grav¬ 
itomagnetic Love numbers of irrotational stars with real¬ 
istic equations of state in Fig. together with the Love 
numbers of strict hydrostatic equilibrium. We note that 
the gravitomagnetic numbers of irrotational bodies tend 
to 0 for realistic equations of state. For completeness, 
we show the gravitomagnetic Love numbers for the poly¬ 
tropic equation of state used in [5] for n = 0.5,1, 2; our 
result agrees perfectly with those of Fig. 1 in [9]. 


III. GRAVITOMAGNETIC LOVE NUMBER 
UNIVERSALITY 
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We have integrated the model equations for 5 represen- 
^ ' tative equations of state, namely APR |53], BSK14 |30j . 
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FIG. 2. The gravitomagnetic Love number for stars in strict 
hydrostatic equilibrium (left) and for irrotational stars (right) 
as a function of the dimensionless moment of inertia (right) 
in a Log-Log plot, for various equations of state. 


AU (called AV14-I-UVII in [21]), SLy [32] and FPS [55] . 
for a range of neutron star mass and radius. We show 
the dependence of the gravitomagnetic Love number for 
strict hydrostatic equilibrium and for irrotational stars 
and stars in strict hydrostatic equilibrium in Fig. [^as a 
function of the dimensionless moment of inertia. 

Our results indicates that both gravitomagnetic Love 
number, for the strict equilibrium and for irrotational 
stars enjoy a quasi universal relations. We find the largest 
deviation between the best fit and our data for irrota¬ 
tional stars (resp. stars in strict equilibrium) is of the 
order of 5% (resp. 10%), and most of the equations of 
state are within « 2% (resp. « 5%). It is in principle 
possible to reduce this error, along the lines of [Mj. We 
did not find a significant enhancement, but we could re¬ 
duce the maximum deviation by almost a factor 2 in the 
case of strict hydrodynamic equilibrium. Note that the 
maximal deviation takes into account less compact stars 
where it is expected that universality deteriorates. In 
all case, we found universality with less than « 5% for 
most of the range of compactness that we have investi¬ 
gated. We use the following form of the function to fit 
the relations 

log{C%) =do+di log(C“/) + d 2 log(C“/)2+ 
d3log(C“/)3 + d4log(C“/)^ (11) 

The coefficients are given in Table [TT| for irrotational stars 
and Table |I] for stars in strict hydrostatic equilibrium. 

Note that the deviation is slightly lower for the irro¬ 
tational case than for the strict hydrostatic equilibrium, 
which is pleasant feature since this case is supposed to 
be more realistic. 

We illustrate the deviations for the irrotational case in 
Fig. I^and for the strict hydrostatic equilibrium in Fig. 
We detail the analysis method we have followed in order 
to try to improve the universal relations in Appendix [A] 



FIG. 3. The dispersion fo the gravitomagnetic Love number 
of irrotational stars The black thick curve is the best fit with 
parameters given in Table [III 



FIG. 4. The dispersion fo the gravitomagnetic Love number of 
stars in strict hydrostatic equilibrium. The black thick curve 
is the best fit with parameters given in Table |I] 


IV. DISCUSSION AND CONCLUSION 

In this short paper, we have investigate the univer¬ 
sality of the gravimagnetic Love number of irrotational 
stars and stars instrict equilibrium, by looking for corre¬ 
lations between the latter and the dimensionless moment 
of inertia, defined by / = I/M^. We have used 5 differ¬ 
ent tabulated equations of state, namely APR, SLy, AU, 
BSK and FPS, that are commonly used in the literature 
on universal relations. 

We found that both Love numbers normalized follow 
approximate universal relation, with a dispersion less 
than « 10% for the strict hydrostatic equilibrium stars 
and less than « 5% for the irrotational stars. Note that 
previous results were reported on the universal behaviour 
of the strict hydrostatic equilibrium case in [25]. We 
found that modifying slightly the normalisation of the 
gravitomagnetic Love number and of the dimensionless 
moment of inertia according to (7“/ and C^k 2 , it is pos- 
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sible to reduce the maximal dispersion by a factor almost 
2 in the case of strict equilibrium. However, the average 
standard deviation did not change significantly. We gave 
the best values of a, b that were determined by minimis¬ 
ing the total standard deviation, obtained by computing 
the deviation to the avarage curve of all the data we 
have generated. We stress that the deviation and mean 
deviation we computed take into account lighter config¬ 
urations, where it is expected that the universality holds 
less. 

We expect this universal relation to extend to a the 
rotation induced quadrupole, as well as to higher grav- 
itomagnetic irrotational multipoles, along with |25j . Fi¬ 
nally, we note that there is a slight improvement in the 
universality when taking into account the internal motion 
of the stars, induced by the gravitomagnetic tidal held. 
We observed than the mean (resp. max) standard devi¬ 
ation in our data is 4% (resp. 14.5%) for the irrotational 
fluid, while it is 7% (resp. 33%) in the case of strict 
static equilibrium. We expect this feature to hold for 
higher multipoles, strengthening the correlations found 
in pS] . 


It should be stressed that the sampling values of 
depend on K. Thus it is not a good idea to use directly 
these datas as they are. Instead, we resampled the data 
on an log-equally spaced sample in the common range of 
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Appendix A: Choosing the normalisation 

Following [34], we investigated a deformation of the 
univeral relation by looking for correlations between 7(7“ 
and k 2 C^, for some a,b to be determined. 

We adopted the following strategy: first we generated 
tables of {IC°‘,k 2 C^) for all the equations of state we 
built for various neutron star configurations. Let us call 
these data and their interpolating function 

where K labels the equations of state. Next we 
have computed the average of the interpolation of these 
data, and obtained a function < labik) >: 

1 ^ 

< Xab{k) >= ^ E (Al) 

K=1 


FIG. 5. logio the cumulative standard deviation of the rela¬ 
tive difference between the curves in the (jfC“, IC^) plan and 
their average, computed with 5 different equations of state for 
irrotational stars. 


the datas. This resampling is now independant of the 
equation of state. Let us call the dispersions computed on 
the resampled data We then compute the variance 
of the joint data, defining 

Dab = Var(U^=ii4f,), (A3) 

where Var is the variance. 

Fig. [^shows \/Dab as a function of a and b. Of course, 
this plot depends on the specific choice of the equations 
of state set, but we expect the minimal region to remain 
roughly the same. As a sanity check, we checked that the 
mean value of the dispersion is a small number, typically 
find of the order of « 10“3. 


where N is the total number of equations of state. 

Then the dispersion of the data of the equation of state 
number K with respect to the avarage values is given by 


/f- < labjKl^) < 
<X„b(Kf) > 


(A2) 


Appendix B: Coefficients of the fit 
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a 

b 

do di d2 ds d^ 

< |A| > Max|A| 

0 

-5 

o o 

-5.07-bO 5.96 - 1 -1.99 - 1 2.58- 2 -1.51 -3 

-8.60-bO 8.69 - 1 -6.53 - 2 2.12 - 3 -2.79 -5 

0.070 0.330 

0.050 0.193 


TABLE I. Values of the parameter of (111, for the direct re¬ 
lation between I and (J 2 and for a different choice of normali¬ 
sation. This case is the case considered in previous literature 
and refered to as strict hydrostatic equilibrium in m- <iAi> 
is the average absolute deviation over all the axis, and Max| A| 
is largest deviation we observe over all the axis. The numbers 
are presented as A ± V = X 10^^. 


a 

b 

do di d2 dz d4, 

< |A| > Max|A| 

0 

-3 

0 

0 

-5.18-bO 1.29 - 1 -7.16 - 2 8.58- 3 -6.47-4 

-5.86-bO 2.80 - 1 -3.38 - 2 1.48- 3 -3.24 -5 

0.040 0.144 

0.034 0.131 


TABLE 11. Values of the parameter of (111, for the direct 


relation between / and k 2 and for a different choice of nor¬ 
malisation. This is the irrotational fluid case considered in 

[9], where internal motion is present. < |A| > is the average 
absolute deviation over all the axis, and Max|A| is largest 
deviation we observe over all the axis. The numbers are pre¬ 
sented as A ± y = A 10^^. 
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